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ABSTRACT 

This paper presents an analysis of core-collapse supernova distributions in isolated 
and interacting host galaxies, paying close attention to the selection effects involved 
in conducting host galaxy supernova studies. When taking into account all of the 
selection effects within our host galaxy sample, we draw the following conclusions: 

(i) Within interacting, or 'disturbed', systems there is a real, and statistically 
significant, increase in the fraction of stripped-envelope supernovae in the central 
regions. A discussion into what may cause this increased fraction, compared to 
the more common type IIP supernovae, and type II supernovae without sub- 
classifications, is presented. Selection effects are shown not to drive this result, and 
so we propose that this study provides direct evidence for a high-mass weighted 
initial mass function within the central regions of disturbed galaxies. 

(ii) Within 'undisturbed' spiral galaxies the radial distribution of type lb and type 
Ic supernovae is statistically very different, with the latter showing a more centrally 
concentrated distribution. This could be driven by metallicity gradients in these 
undisturbed galaxies, or radial variations in other properties (binarity or stellar 
rotation) driving envelope loss in progenitor stars. This result is not found in 'dis- 
turbed' systems, where the distributions of type lb and Ic supernovae are consistent. 

Key words: Galaxies : interactions galaxies : fundamental parameters - galaxies : 
starburst supernovae : general 



1 INTRODUCTION 

The host galaxies of core-collapse supernovae (CCSNe) give 
crucial information when trying to determine the progen- 
itors of these explosions, and have led to some important 
breakthroughs in the study of CCSNe. These include the al- 
most exclusive presence of CCSNe in l ate-type galaxies (e.g. 
Ivan den Bergh. Li fc Filippenkoll2005r ) implying progenitors 
which only occur in young stellar populations. More recently, 
host properties played an important role in the discovery of 
a new classification of supernovae termed 'calcium-rich', and 
showing spectroscopic characteristics of CCSNe yet origina t- 
ing in regions devoid of star-formation (jPerets et al.| [2010). 

Even more information can be gathered on the possi- 
ble supernova progenitors of the various subtypes of CC- 
SNe, by looking at the specific site of the SN explosion 
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within the host galaxy. The CCSNe group divides broadly 
into type II (SNII) and types lb and Ic (SNIbc), dis- 
tinguished by the p resence or lack o f hydrogen in their 
spectra respectively (Minkowski 1941). The SNII subclass 
can be further divided into IIP, IIL, Iln, lib which will 
be discus sed later; f o r a r eview of supernova classifica- 
tions sec Filippcnk^ l| 19971 ). Many studies have investi- 
gated the association of these CCSNe subtypes with spiral 
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These analyses indicated that SNIbc are both more strongly 
associated with HII regions implying higher mass progeni- 
tors, and more centrally concentrated within galaxies, when 
compared to type II SNe, indicating a metallicity depen- 
dence with observed SN-type. Both interpretations can be 
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explained theoretically by mass-loss via radiatively-driven 
stellar winds whose effe ctiveness increases with progenitor 
mass (IHeeer et al.ll2003ri. or high meta llicity environments 
Ce.g. lPuls eta l. 1996; H eger et al.l l2003h- 



In 
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HAJ10, the current authors carried out a study on the dis- 
tributions of CCSNe types in both 'undisturbed' and 'dis- 
turbed' host galaxies. This study concluded that within 
the disturbed systems an excess of SNIbc existed within 
the central regions when compared to those of type II. 
A case study on Arp 299, a so-called 'supernova-fa ct ory ' 
i|Neff. Ulvestad fc Tendl2004:lPerez-Torres et al.ll2009h, also 



carrie d out by the authors IjAnderson. Habergham. fc James! 
l201ll henceforth All), found the same unusual central ex- 
cess of stripped-envelope supernovae (SE-SNe). These re- 
sults have ignited a discussion into the role selection effects 
can play within host galaxy studies. 

The analysis within this paper will concentrate on SE- 
SNe against SNII+IIP (consistent with All); where SNII 
are all type II supernovae without sub-classifications, and 
SE-SNe include all type lb, Ic, Ib/c, lib and ILL. Any type 
Iln supernovae (SNIIn) within the sample have been omit- 
ted from the analysis and discussion due to uncertainties in 
their progenitor stars , and often in their classification (e.g. 
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The accumulation of subtypes into the SE-SNe group, 
rather than the SNIbc and SNII groups of HAJ10, comes 
as a result of recent studies into the progenitors of CCSNe 
and the e mergence of a sequ ence of progenit or envelope- 
stripping (iDessart et al.l l201l[). or mas s loss (| Heger et al.l 
bOOilCrowtherlteOOTHGeorgv et alj|2009h . The sequence has 



type IIP retaining their entire stellar envelope, followed by 
an increasing loss of stellar envelope material from IIL, to 
lib, lb and Ic, which have lost both their hydrogen and 
helium layers. The interpretation for this increased enve- 
lope stripping has been debated. Theoretical studies suggest 
that it could be due t o an increase in the progenitor mass 
fe.g. iHeger et al.ll2003r ). the metallicity of the environment 
fe.g. lHeger et a l. 2003), t he presence of a binary companion 
(e.g. IPodsiadlowski. Joss, fc Hsul Il992h . or stellar rotation 
(e.g. iMevnet fc Maederj I2003TT Observational evidence for 
any one particular theory is difficult to determine, though 



there is evidence for both metallicit y (Anderson & James 



2009 hereafter AJ09, iModiaz et ail l201ll 
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201lf ). and progenitor mass ( Anderson fc Ja mes 2008) play- 



ing a substantial role in the stripping of stellar envelopes 
prior to explosion. The best observational constraints would 
come from direct detections of CCSNe progenitors on pre- 
explosion images. Although this te chnique has had great 
success in a number of cases ( e.g. [p ins- Rosa et al.l l201ll; 
iMaund et al] 120111 . see ISmarttj 120091 for a review on the 
topic), the need for very nearby events severely limits the 
statistics available from these studies, which thus far have 
only had success with the much more common type II sub- 
classes. In reality the cause is likely to be a combination of 
all of these effects. 

The nature of all CCSNe, having high-mass, short-lived 
progenitors, provides a direct probe of recent or on-going 
star-formation within their host galaxies. The differences 
between the different subtypes also allows constraints to be 
placed on this star-formation. So investigations into the host 



galaxies of CCSNe provide insights both ways: into the pos- 
sible progenitors of the supernovae themselves, and, into the 

mode of star formation in the host. 

Our previous papers, culminating in lAnderson et al.l 

(2012), have established a probable mass sequence for CC- 
SNe progenitors. Building on this the aim of this paper 
is to test the robustness of the authors' previous study 
into the variation of CCSNe distributions with environment 
(HA J 10). We will present the data and host galaxy study 
sample in Section 2, followed by a new analysis of the ra- 
dial distribution of CCSNe within their hosts in Section 3. A 
rigorous investigation into possible selection effects present 
within the sample will be discussed in Section 4. Section 5 
will explore possible interpretations of the results, namely 
the effects of SNe environment metallicity, and the overall 
metallicity gradients of the hosts; the contribution of stellar 
rotation and binarity to the possible progenitor population; 
and the possibility of a modified initial mass function (IMF) 
preferentially producing the most massive stars in the cen- 
tral regions of 'disturbed' host galaxies. Finally, in Section 
6 we will summarise the results and present our preferred 
interpretation and any possible implications. 



2 DATA SAMPLE 

Our data are drawn from observations of host galaxies indi- 
cated by SN discoveries contained within the Padova-Asiago 
SN catalogue^] and IAL0 SN catalogues. Such catalogues in- 
herently contain biases with some resulting from initial SNe 
studies which targeted bright star-forming galaxies, rather 
than the 'blind' searches many groups are cur rently pur- 
suing (e.g. Palomar T ransient Factor y (PTF), iRau et al.l 
120091 ; lLaw et al.l 120091 ; Pan-STARRS, iKaiser et all l201Ch . 
Our analysis contains 218 host galaxies containing 280 CC- 
SNe of all subtypes, and although biases still persist within 
this sample, until such sufficiently large samples are accu- 
mulated through 'blind' searches, this is the best analysis 
possible. The large sample within this study also enables us 
to cover more of the rarer subtypes of SNe, but given the 
selection effects contained in the sample we do not claim 
that the relative numbers found here are representative of 
the Universe as a whole, and the conclusions drawn from the 
data are not dependent upon statistical completeness. 

These host galaxies have been observed over 15 years 
of observations in both SNe-related and Ha galaxy studies. 
This provides us with a heterogeneous and randomly sam- 
pled selection of the to-date, observed local (median host 
recession velocity ~1918 kms -1 ) CCSNe. 

All galaxies have Ha narrow-band and 7?-broad band 
imaging, and come from the following facilities: the 
MPG/ESO 2.2m telescope at La Silla, Chile; and the Liver- 
pool Telescope (LT), Jacobus Kapteyn Telescope (JKT) and 
Isaac Newton Telescope (INT) all on La Palma, the Canary 
Islands. The data were reduced in a standard manner using 
routines in Starlink and IRAF0. For more detailed informa- 
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3 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
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tion on the observations and data reduction process refer to 
AJ08 and Anderson et al. 2012 (submitted). 

When analysing host galaxies by a characteristic such 
as disturbance, classification errors can become problematic. 
In HAJ10 the classification of galaxy disturbance was car- 
ried out by two independent visual inspections of the host 
galaxy images. Although more quantitative methods exist 
for classifying galax y disturbance, namely through the de- 
gree of asymmetry (jConselice. Bershadv. fc Jangrer] l2000l ; 
lLotz. Primack fc Madau" 20041 ) . the classification within 
this study assesses such a wide range of parameters that 
it is better done via an eye-ball study than any one quan- 
titative singl e measu r e. Th e vis ual classification schemes 
descri bed in ISuracd (jl998r ) and IVeilleux. Kim fc Sanders 
(120021) have been widely us ed (e.g. iMiralles-Caballero et~ 
120111 ; IZamoiski et aljkulll ). and in this paper we adopt a 
similar technique, attributing galaxy characteristics to lev- 
els of disturbance. Although the most disturbed galaxy sys- 
tems often also show signs of merger-triggered starbursts 
in the nuclei, there are undoubtedly some galaxies which 
show signs of interaction without displaying such energetic 
star-formation. For the reanalysis presented in this paper 
the authors themselves reclassified each galaxy within the 
sample as 'disturbed' or 'undisturbed' based on a number 
of characteristics. An 'extreme' sample of disturbed galaxies 
was chosen, which contained hosts undergoing a major in- 
teraction with another galaxy, or a which had large degree 
of asymmetry indicating a recent tidal interaction with an- 
other system. A larger sample of 'disturbed' galaxies (includ- 
ing the 'extreme' systems above) included those displaying 
at least two signs of minor interactions. These signs were; 
presence in a group with minor interaction, minor degrees 
of asymmetry, shells, tails, irregular structure (small irregu- 
lar systems are not included in this analysis - this irregular 
structure is within a larger star- forming galaxy), a double 
nucleus, extreme or irregular dust patterns or a clumpy mor- 
phology. 

A full list of classifications for each galaxy in our sample 
is presented in Table 1. The full table is available online 
through the supplementary material. Q 
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Figure 1. Fractional i?-band light distribution of CCSNe in 
'undisturbed' host galaxies. 

• Undisturbed: 0.766 (59 : 77) 

• Disturbed: 1.14 (65 : 57) 

• Extreme: 1.85 (26 : 14) 

It is clear that the disturbed galaxies have a higher frac- 
tion of SE-SNe, by more than a factor of two, when com- 
pared to undisturbed galaxies. 

We now present the analysis of the radial distribution 
of 258 CCSNe (the 22 type Iln SNe within this sample will 
be the subject of a future paper), in both 'undisturbed' 
(136 SNe in 114 hosts) and 'disturbed' hosts (123 SNe in 
89 hosts). The radial distribution is defined in terms of the 
fraction of _R-band light Fr(iJ) which lies within the circle 
or ellipse which contains the SN. This means that an Fr(R) 
value of 0.0 corresponds to a supernova at the central R- 
band peak of the galaxy emission, while a value of 1.0 implies 
an extreme outlying SN. A full description of this analysis 
is given in AJ09. 



3 RESULTS 

We begin this section by testing whether the ratio of SE-SNe 
to SNII differs with galaxy disturbance. The absolute rates 
of each SNe sub-type would be subject to bias within our 
sample regarding host galaxy selection. However, between 
the different galaxy samples no such bias is present and so 
the relative frequencies of each CCSNe type give us useful 
information. Below we list the ratio of SE-SNe to SNII 
in each sample, with the numbers of SNe in each group 
presented in brackets. 

SE-SNe : SNII+IIP 

for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 

4 Where starbursts have been indicated these come from the pres - 
ence of the galaxy within the GOALS sample l lArmus et al. 2009) 
as either a LIRG or ULIRG, or where the system is noted as a 
starburst on the NASA/IPAC Extragalactic Database (NED). 



3.1 Radial Distribution Analysis 

In Figure 1 we present the CCSNe distribution in the iso- 
lated or 'undisturbed' galaxy sample, and in Fig. 2 in the 
'disturbed' sample. The radial distributions are presented 
as histograms of the Fi(R) value in which the top panel 
presents all CCSNe, SNII (type IIP and SNII without a sub- 
classification) in the middle panel and all SE-SNe (lb, Ic, Ibc, 
IIL and lib) in the bottom panel. The combination of SNIIP 
with SNII without sub-classifications is based on the as- 
sumption that ~70% of these SNe would be SNIIP based on 
Lick Observatory Su pernova Search (LOSS) observed frac- 
tions (|Li et al.l l201ll L Indeed a Kolmogorov-Smirnov (KS) 
test of SNII+SNIIP against SNIIP within the whole sam- 
ple gives a probability that they are drawn from the same 
parent distribution of ~86%. 

As discussed in Section 2, the classification criteria used 
to define these samples mean that within the 'disturbed' 
sample are galaxies in all stages of interaction. This may 
be a galaxy within a group with signs of minor asymmetry, 
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Table 1. An excerpt of the galaxy classification information for our sample. The full table is available online. 
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Figure 2. Fractional ij-band light distribution of CCSNe in 'dis- 
turbed' host galaxies 

to merging galaxy pairs. A subset of the 'extreme' galaxy 
sample can be seen in Fig. 3, and the results of the radial 
analysis of the CCSNe in these systems in Fig. 4 with the 
same panels as Figs. 1 and 2. 

The distributions of SNII+IIP and SE-SNe within each 
of the samples appear markedly different. We present tests 
on the statistical significance of these differences using a 
KS test. The KS test takes two parameters to calculate the 
probability; the 'distance' (the largest difference between 
the distributions on a cumulative distribution plot) between 
the distributions, and the number of events within each 
distribution. Due to the decreased number of events in 
the 'extreme' sub-sample throughout the analysis of these 
results we will present both the probability (P) and the 
distance (D) of the distributions, in order to establish 




Figure 3. Mosaic of 12 of the 28 host galaxies contained within 
our 'extreme' sample. 

whether a reduced significance in 'P' is due to either a 
smaller sample size or a change in the distribution. 

The results of the KS tests carried out in this study are 
presented in Table 2, and will be discussed throughout the 
rest of this section. 

Within the 'disturbed' galaxy sample there exists a statisti- 
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Figure 4. Fractional R-band light distribution of CCSNe in the 
'extreme' host galaxy sub-sample 
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Figure 5. Fractional /J-band light distribution of type IIP and 
Ibc SNc in 'undisturbed' hosts. 



cally significant (P~6 per cent) difference between the dis- 
tributions of SE- and IIP- supernovae; a result which be- 
comes more exaggerated within the 'extreme' sample. This 
is a result that stands from HAJ10 despite the addition of 
all SE-SNe into the same class. 

In Figs. 5 and 6 we present the same data as in Figs. 
1 and 2, but without the intermediate SE-SNII classes 
(IIL and lib), and without type II SNe without sub- 
classifications, for 'undisturbed' and 'disturbed' galaxies re- 
spectively. In both figures the top panel shows the overall 
distribution of all CCSNe in the systems, the middle panel 
the SNIIP and the bottom panel the SNIbc. The distribu- 
tions shown therefore represent the most extreme within our 
sample, where only SNe which have lost their entire hydro- 
gen envelope are plotted, and the uncertainties in the SNII 
unclassified sample are removed. 

Although the statistics are clearly reduced in these 
examples, the central excess of SNIbc becomes even more 
apparent in the 'disturbed' galaxies. Despite the small 
statistics a Kolmogorov-Smirnov test places the probability 
that SNIIP and SNIbc are drawn from the same parent 
population within disturbed systems at ~4 per cent. The 
results of the KS tests are shown in Table 2. 

Combining the statistically consistent SNII popula- 
tion with the SNIIP the results become more dramatic. 
Within the 'disturbed' galaxies there is only a ~1 per cent 
probability that the SNII+SNIIP and SNIbc populations 
have the same parent source. The reduced statistics in the 
'extreme' sample raises this probability to ~3.5 per cent, 
although the large 'D' value indicates that the distributions 
of these populations are significantly different (Table 2) . 
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Figure 6. Fractional i?-band light distribution of type IIP and 
Ibc SNe in 'disturbed' hosts. 

3.2 Distribution of SNIb and SNIc in undisturbed 
galaxies 

Where the statistics allowed we explored the distributions of 
all CCSNe types, and in doing so another significant result 
has emerged. Our samples are now large enough so that 
the SNIbc class can be broken down into SNIb and SNIc 
individually. Within the 'undisturbed' sample we have 21 
SNIb and 23 SNIc, and in the 'disturbed' sample, 21 SNIb 
and 27 SNIc. Figures 7 and 8 display the histograms of the 
distributions of each population, where the top panel is again 
the overall distribution of CCSNe, the middle panel SNIb 
and the lower panel SNIc. 

The principal result here is that following a KS test, 
although the distributions of SNIbc within the 'disturbed' 
sample are completely consistent (P~20 per cent), within 
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Table 2. Kolmogorov-Smirnov test results for Section 3. 
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the 'undisturbed' sample the two distributions are very 
different. 

This places a probability that SNIb and SNIc, within 
'normal' spiral galaxies, are drawn from the same parent 
population of only ~0.4 per cent (Table 2). Although the 
statistics here are small this is an astonishing result, and 
not one which can be easily explained in terms of selection 
effects. 

The result seems to be driven in part by the lack of 
SNIb in the central regions of these systems, and in part 
by the lack of SNIc in the outer regions of 'undisturbed' 
galaxies. The probability of 'missing' SNIc explosions in the 
outer regions of galaxies, where the host extinction is low, 
seems unlikely. Correspondingly , the similar p eak absolute 
magnitudes of SNIb and SNIc ()Li et al.l[201ll ). mean that 
it is unlikely SNIb are being missed in the central regions, 
where SNIc are still observed. 

The most likely explanation for the different distribu- 
tions of SNIb and SNIc in 'normal' host galaxies is that the 
metallicity gradient within the host plays a major role in 
evolving the progenitor star prior to explosion. Envelope- 
stripping through radiatively driven, metallicity-dependent, 
winds has been discussed in Section 1. We believe the re- 
sults presented here, showing that the distributions of SNIb 
and SNIc within 'undisturbed' galaxies are statistically very 
different, provides direct evidence for this. 

One interesting aspect of our results is the apparent lack 
of SNIc in the outer regions of our 'disturbed' galaxies, where 
we would expect som e of the lowest meta llicity regions. This 
echoes the findings of lArcavi et al.l (l2010h who find a lack of 
SNIc in intrinsically low metallicity dwarf galaxies. 

The contribution of metallicity within 'disturbed' host 
galaxies does not appear to have as strong an effect as we 
will demonstrate in Section 5.1. 



4 SELECTION EFFECTS 

Within any host galaxy supernova study, selection effects 
and biases can enter data from two possible sources: the 
supernova selection, and the host galaxy selection. We shall 
explore the different effects arising from each in this section. 

4.1 Supernovae 

A range of selection effects can contribute to biases within 
a CCSNe study. The first of these is that different CCSNe 
types have different magnitude ranges, both in peak magni- 
tude and the duration of the light curve, and so within any 
study, one is likely to lose intrinsically fainter SNe subtypes, 
which are never discovered due to the brightness of their host 
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Figure 7. Fractional i?-band light distribution of lb and Ic SNe 
in 'undisturbed' host galaxies. 
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Figure 8. Fractional ij-band light distribution of lb and Ic SNe 
in 'disturbed' host galaxies 

galaxy, or the extinction within it. It is often accepted that 
SNIIP are most likely to be unaccount e d for in large volum e 
galaxy studies (|Pastorello et all |2004 ISmartt et ail 120091 ). 
due to their fainter peak magnitudes during outburst. This 
effect could be exaggerated within 'disturbed' host systems 
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due to the increased dust content, and warped morphology. 
Howev er, careful ana lysis of the SN luminosity function from 
LOSS (|Li et al.ll201lh . found no indication that SNIIP have 
intrinsically fainter explosions. Although SNIIP have a faint 
tail (not unlike SNIc) , the absolute magnitudes of each CC- 
SNe subclass have enough overlap that missing SNIIP are 
unlikely to explain the results pre sented in this paper. 

The Shaw Effect (|ShawHl979h , that it is more difficult 
to detect supernovae in the inner regions of distant galaxies, 
also affects all supernovae searches and hence all host galaxy 
sample analysis. Although the Shaw effect was much more 
applicable to photographic plate searches, where the centres 
of galaxies were often over-exposed, it is still a source of bias 
in today's SN searches. Studies have been carried out into 
the number of SNe estimated to be lost through the Shaw ef- 
fect. ICappellaro et afl l|l993h found that out to host galaxy 
recession velocities of ~6000 kms -1 approximately 35 per 
cent of all CCSNe types may be lost in photographic plate 
searches in the inner regions of host galaxies. This compares 
to ~22 per cen t of all SNe out to ^4000 kms -1 in som e of 
the first visual (Evans, van den Bergh fc McClur ell 19891 ) and 
CCD based automated searches (jMuller et al.lll992l ). Within 
this range of host galaxy redshifts there is no preferential 
bias to any particular type of SNe though statistics at this 
time were poor. Studies into the Shaw effect on modern SNe 
samples, assembled using modern CCDs and both targeted 
and 'blind' searches are lacking, and so f or the purpose of 
this in vestigation the results drawn from ICappellaro et afl 
|l993t ) are the best estimates of bias available in the litera- 
ture. 

To test further for selection biases, we analysed 
the distribution of discovery magn itudes (taken from 
iLennarz. Altmann &i Wiebuschl |2012| and the Padova- 
Asiago SN catalogue 0) for our sample of SNIIP+II and 
SE-SNe. All SNe with available discovery magnitudes were 
colour-corrected to the it-band, due to this being closest to 
unfiltered light (a large proportion of the discovery magni- 
tudes). Our corrections were based upon available CCSNe- 
subtype multi-band photometry on the SUSPECT database 
(the online Supernova Spectrum Archive). 

For this analysis we used the discovery magnitude of 
each SN rather than a peak magnitude. In the context of 
exploring selection effects within our sample, specifically 
whether we are missing 'faint' SNe, the discovery magnitude 
is the most important. Figure 9 shows the absolute _R-band 
magnitude at discovery of every CCSNe within our sample, 
plotted against redshift. Here SE-SNe are represented by 
the blue open circles, and SNII+IIP the red crosses. As ex- 
pected, there is a much larger range of discovery magnitudes 
within our sample at low redshift. However, it is interesting 
to note that the brightest events are found at all redshifts, 
and not just in the most distant galaxies. 

The importance of analysing the discovery magnitude 
of each SNe against the Fr(R) value used in this paper is 
paramount in establishing any possible selection effects. If 
we were to lose 'faint' SNe, we would do so in the central, 
dusty and high surface brightness regions, defined in this pa- 
per as within ^20 per cent of the _R-band light (Fr(i?)<0.2). 
Figures 10 and 11 display the distribution of absolute and 
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Figure 10. Fractional _R-band light distribution of each CCSN 
against its absolute magnitude at discovery. 



apparent discovery magnitudes, respectively, of each SN in 
our sample against their Fr(i?) value. Again, SE-SNe are rep- 
resented by blue, open circles, and SNII+IIP by red crosses. 

Both Figs. 10 and 11 show no correlation between the 
discovery magnitude (absolute or apparent) and the Fr(i?) 
value, for either the SE-SNe or the SNII+IIP. Within the 
central regions we observe all CCSNe types, over a similar 
range of magnitudes. 

Further to this, it is important to analyse whether all of 
the central SNe within our sample are present at the lowest 
redshifts. It may be possible that the central SNe we detect 
are all in the most local galaxies, and hence a selection ef- 
fect could still affect the result. In order to analyse this we 
have plotted the recession velocity of each host galaxy and 
the Fr(R) values of the SNe within it. This is displayed in 
Fig. 12, where again SE-SNe are the blue open circles and 
SNII+IIP are the red crosses. The dashed black lines repre- 
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Figure 11. Fractional i2-band light distribution of each CCSN 
against its apparent magnitude at discovery. 



sent the median recession velocity of our sample and 6000 
kms , which was the limit of the HAJ10 study, and the 
limit to which we would claim our results are reliable, due 
to the small number statistics above this range. 

Figure 12 again shows that there is no correlation 
between the Fr(_R) value of the SNe and its host galaxy 
recession velocity. Certainly out to ~6000 kms -1 there is no 
indication of losing any II+IIP within the central regions, 
or at higher distances from the centre. A KS test between 
low and middle redshift bins shows that the distribution of 
SNII+IIP are completely consistent (P=0.990, D=0.0814). 
The overall distributions of CCSNe within the 3 redshift 
bins (z<median; median>z>6000 kms -1 ; and z>6000 
kms ) are all consistent with being drawn from the same 
parent population. 

We next test the radial distributions of CCSNe of all 
types within the different redshift bins. 

KS test: 

• low z bin vs middle z bin: P=0.947 D=0.0647 

• middle z bin vs high z bin: P=0.193 D=0.2353 

This shows that these Fr(R) distributions are statisti- 
cally completely consistent and that the main result pre- 
sented in this paper, of an excess of SE-SNe within the cen- 
tral regions of 'disturbed' hosts, is not due to a failure to 
detect SNIIP within these regions. 

Despite the small number of events in our sample above 
6000 kms -1 , all of the histograms provided in Section 3, 
and all of the KS tests given, include these events. Although 
the sample at this redshift range is far from complete, the 
inclusion of such objects within the results have no effect on 
the outcome, as the distribution of events at all redshifts is 
essentially flat and therefore cannot drive the results found. 

In a very conservative test of possible distance depen- 
dent selection effects, we now minimise these by analysing 
only those within 2000 kms -1 . This reduces the sample size 
in both 'undisturbed' (67 CCSNe) and 'disturbed' (56 CC- 
SNe) samples. Even in this sample, a central excess of SE- 
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Figure 12. Fractional _R-band light distribution of each CCSN 
against the host galaxy recession velocity. The lower dashed line 
represents the median recession velocity for the sample, and the 
upper one the 6000km- 1 cut. 



SNe is still visible, providing 11 of the 17 events within 
Fr(i?)=0.2 (compared to 3 out of 7 in the 'undisturbed' sys- 
tems). However, at this level the result is no longer statisti- 
cally significant. 

Our radial analysis can also be carried out on the Ha 
emission within a galaxy, which traces star formation rather 
than the old stellar light (traced by the -R-band). Here an 
Fr(ifa) value of 0.0 means the supernova is closer to the 
central .R-band peak of the galaxy than any Ha emission, 
and again a value of 1.0 indicates an outlying SN. 

Figure 13 shows the histogram of Fr(ifa) in our 'dis- 
turbed' galaxy sample. This figure shows that the central 
distribution of SE-SNe is apparent, even with respect to the 
Ha emission within the host galaxy. This suggests that it is 
not due to a selection effect. 

All of the investigations carried out in this section sug- 
gest that the results found in Section 3 cannot be driven by 
any selection effect. We find in our data an absolute excess 
of stripped-envelope SNe within the central regions of 'dis- 
turbed' galaxies, even with respect to the star formation as 
traced by Ha. 



4.2 Host Galaxy Sample 

It is well documented that galax y disturbance is linked 



it is well documented that galaxy disturbance is luiKcd 
to rapid bursts of star formation (lLarson fc Tinslev 19781: 
Joseph et aDll984lKennicutt fc Keellll984lKennicutt et all 



1987h which are frequently concentrated in the central re- 



gions o f the distur b ed or merging systems dJoseph fc Wright] 
Il985l : iKeel et all Il985l ). probably linked to the cen- 
tral gas concentrations seen in merger simu lations 
l|Barnes fc Hernquisdll99ll : lMihos fc HernauistJll996l ). Thus, 
a possible reason for the excess of central stripped-envelope 
SNe in disturbed galaxies is that these regions are domi- 
nated by starbursts of such extreme youth that they have 
not yet started to produce SNIIP (under the assumption 
that CCSN type is driven primarily by progenitor masses, 
and hence correlate strongly with stellar lifetimes). 
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Figure 13. Fractional Ha light distribution of CCSNe in 'dis- 
turbed' hosts. 

In extremely youn g star formation regions, suc h as some 
observed in Arp 299 (jAlonso-Herrero et al.l l2000h , it is to 
be expected that current supernovae should have high mass 
progenitors, regardless of the form of the IMF. However, this 
requires the galaxies to have been caught at just the right age 
after an extremely short epoch of starburst activity, which 
requires some fine tuning. 

The degree of fine tuning required depends on the 
ranges of initial mass of progenitors resulting in different 
types of supernovae, and in particular the initial mass corre- 
sponding to the transition between SNIIP and the stripped- 
envelope types. This is a very uncertain value, and depends 
on other factors, e.g. metallicity and rotation. Estimates 
for the transitional mass between a type II and SE-SNe 
for s i ngle-star progenitors i nclude 25-40 M m jHeger et al" 
20031 ; lEldridee fe Toutll20oi ) and ~ 25 M (|Gal-Yam et al 
2007 j ). However, th e rotat ing st ar models pres e nted r ecently 



by lEkstrom et all (|2012T ) and iGeorgy et all (|2012h lower 
this value to between 15 and 20 Mq. Observationally, the 
highest masses attributed to SNIIP progenitors that have 
been detec ted in pre-explosion im aging are 15-18 M for 
SN1999ev ilMaund fc Smarttl 120051 ) ; 15 M for SN2004dj 



Mai'z-Apellaniz et al.ll2004 ); and 15 M again for SN2004et 
Li et al.l l2005i r although a recen t study of the progenitor 



of SN2012aw (|Fraser et alj|2012h . an unconfirmed SNIIP, 
found a higher mass range of 14-26 M . From the statistics 
of all s uch progenitor detections found at the time of publi- 
cation, ISmarttl (2009) derived an upper limit to the masses 
of SNIIP progenitors of 16.5 M . Thus there is strong evi- 
dence that stars at least as massive as 15 M can produce 
SNIIP, and conservatively taking this as the maximum mass, 
a burst of star formation would start producing SNIIP after 
some 13 Myr, accord ing to the high-mass stellar lifetimes of 
iMevnet et al.l (l994|). If instead the theoretically-preferred 
higher mass limits for this transition are adopted, SNIIP 
could be occurring in as little as 5 Myr. While it is possible 
to catch one galaxy in a sufficiently recent starburst phase 
that no SNIIP are yet occurring, for an ensemble of galaxies, 
such as those presented here, and previously in HAJ10, it 
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Figure 14. Examples of the 'disturbed' supernova-hosting galax- 
ies from this study, arranged by the inferred stage of the 
merger/interaction. Merger timescale increases from ~260 Myr 
at top-left to ~910 Myr at bottom-right. 



seems extremely improbable that most could be caught at 
such a critical phase of activity purely by chance. 

One possibility is that our selection by morphological 
disturbance may have resulted in systems with the very 
young starbursts required for this timescale explanation. In 
an attempt to determ ine the plausibility o f this explan ation, 
we used the results o f lMihos fc Hernauistl (|l994Ul996h . who 
used Smoothed Particle Hydrodynamics simulations to fol- 
low the behaviour of gas within minor and major mergers 
of galaxies. These models contain a prescription for star for- 
mation based on t he original volu me density formulation of 
the Schmidt Law (|Schmidtlll959l ). enabling the calculation 
of the star formation rate as a function of time. The over- 
all timescale is defined by a system of units that assumes 
the original galaxies, prior to merging, to have had proper- 
ties similar to the Milky Way. This is a reasonable approx- 
imation for the bright interacting and merging systems of 
interest in the present paper. 

The details on the classification of galaxies within this 
study were presented in Section 2 and Table 1. Figure 14 
shows a subset of those galaxies classified as 'disturbed', 
which also hosted centrally-located CCSNe. These have been 
matched up by eye wit h the stellar mass distribu tions of 
the simulated mergers in Mih os fc Hernauistl l|l996l ) and ar- 
ranged in an approximate chronologic al order, according 
to the corresponding time within the iMihos fc Hernauistl 
jl996l ) simulations. For the 6 systems shown, the resulting 
merger ages are 260, 260, 689, 780, 845 and 910 Myr re- 
sp ectively, from top l e ft to bottom right. The simulations 
of lMihos fc Hernauistl (|l996l ) also reveal the stages at which 
the central gas densities, and hence the star formation activ- 
ity, are predicted to be strongly enhanced. This is the case 
for all but the first two frames in Fig. 14, with the predicted 
starburst phase covering a simulated time interval of over 
200 Myr. 

While the exact ages derived from this comparison of 
merger morphologies with simulations are subject to sub- 
stantial uncertainties on a case by case basis, the overall 
timescales implied by the differences in galaxy morphology 
for the observed sample cover a range of hundreds of Myr, 
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far in excess of the finely-tuned ages of a few Myr required to 
select the 'stripped-envelope dominated' phase at the start 
of a starburst. Even in the idealised case where a merger- 
induced starburst can be represented by a narrow spike in 
the central star formation rate, then for every one merger in 
the narrow age window where only stripped-envelope SNe 
are being produced in these regions, there should be several 
somewhat older starbursts producing only SNII. There is no 
observational evidence for the latter phenomenon. 

Apart from this fine-tuning argument based on the 
appearances of galaxies and their dynamical timescales, 
there is further independent evidence that observed star- 
bursts have maintained their activity for timescales easily 
long enough to be producing all types of CCSNe. Stud- 
ies of the energetics and spectroscopic properties of star- 
bursts generally infer ages of ^40 Myr for the s tarburst 
phase jSanders fc Mirabellll996l ; ICaputi et alj|2009h . Of the 
most direct relevance for the present argument is the detec- 
tion of extended sof t X-ray emission from the regions sur- 
round ing starbursts (|Franceschini et al]|2003l ; llwasawa et al.l 
This emissio n appears to be driven by CCSNe 
l lwasawa et aljfeoilf ). and the extent of the emission im- 
plies that these SNe have been occurring for ~10 8 years 
(K. Iwasawa, private communication). Again this is easily 
long enough to have reached an equilibrium state in which 
all types of CCSNe should be fairly represented. 

Therefore we conclude here that the selection of young 
starbursts does not drive the result presented here or in 
HAJ10. 



5 INTERPRETATIONS 

Given that we have explored all possible sources of selection 
effect and established that none can drive the central excess 
of SE-SNe in disturbed galaxies, we now explore the possi- 
ble interpretations of this result. We look at each possible 
channel of progenitor envelope-stripping, and how the en- 
vironment within the central regions of interacting galaxies 
may affect each channel. 



5.1 Metallicity 

As discussed briefly in the introduction, one key parame- 
ter thought to lead to envelope-stripping, and hence the 
observation of a SNII or SNIbc, is metallicity. Many SN 
host studies have observe d a central concentration of SNIbc 
when compared to SNII (Bartunov, Makarova, & Ts vetkovl 
19921 ; Ivan den Berglj[l99~7i; iTsvetkov. Pavlvuk fc Bartunovl 



2004; lHakobvan et al.1 120091 ; iBoissier fc Prantzoj 1200a 
Anderson fc James! 12009] ). as seen here, and this has gen- 



erally been assumed to be the result of metallicity gradi- 
ents present within the host galaxy. Within a 'normal' star 
forming galaxy, metallicity gradients have been observed, 
with the central re gions being more metal rich than the 
outer regions (e.g. iHenrv fc Worthevl 1 19991 ). Within this 
scenario central SNe are more likely to lose their outer 
envelopes through metallicity-dependent line driven winds 
JPuls et al l Il99d : iKudritzki fc Puis! |2000| ; iMokiem et al.l 
12007ft . All SNe host galaxy studies to date have failed 
to distinguish between disturbed or interacting systems 
and 'normal' spirals within the sample (although see 



iPetrosian et al.ll2005l for a study of CCSNe with relation to 
starbursts). Within merging galaxies any metallicity gradi- 
ent originally present is disrupted by the interaction and 
is smoothed out within the system dKewlev et al.l 120101 ; 
iRupke. Kewlev fc Chienl l20ld ; iRich et al.1 |2012| ). through 



the infall o f pristine gas from the outer regions into the 
centre (e.g. Rampazzo et al.l 20051; Hibbard fc van Gorkoml 



ll99fj ; lKewlev, Geller fc BartonlbOOet). often fuelling a burst 
of star formation (|Barnes fc Hernq uist 1996). Furthermore, 
studies of the metallicity of SNe sites have found little dif- 
ference between regions that host SNII and SNIbc, and a 
presence of both majo r groups of SNe at all metallicities 



Anderson et al. 2010h. This is supported by a study of 



NGC 2770 (|Thone et ail 120091 ') which found that the large 
number of SNIb within the host were found in low metallic- 
ity regions, negating envelope stripping via radiative winds. 
These studies suggest that the central excess of SE-SNe in 
disturbed galaxies is not due to metallicity. 

In order to study any metallicity gradient present within 
this sample, we obtained spectroscopy for a random subset 
of the host galaxies. Long-slit spectra for 23 of the host 
galaxies were obtained using the Intermediate Dispersion 
Spectrograph (IDS) on the INT on La Palma, in the Ca- 
nary Islands. The slit was positioned to pass through the 
nucleus of each galaxy, to give a positional reference, and an- 
gles chosen to coincide with as many HII regions as possible. 
Met allicities were determine d using the line ratio diagnostics 
from IPettini fc Pagell l|2004 ) , who use emission lines close to- 
gether in wavelength to negate extinction effects, which also 
addre sses any atmospheric differential refraction (Filippcnko 
Il982l ) issues. In any case, observations were taken where pos- 
sible at the parallactic angle. It was not always possible to 
obtain a spectrum at the SN-host HII region doing this and 
so where possible (at low airmass) a second angle was ob- 
served to obtain metallicity measurements for specific SNe 
sites. Standard reduction procedures were employed using 
Starlink packages, and spectra extracted from each observ- 
able HII region along the slit. Each spectrum was then wave- 
length and flux calibrated. A subset of the derived metallic- 
ity gradients is presented in Fig. 15. Each column in Fig. 15 
represents a different sample within the current study; the 
first column displays host galaxies in the 'undisturbed' sam- 
ple, the middle column 'disturbed' galaxies, and the final 
column galaxies within our 'extreme' subset. Any galaxies 
containing known Active Galactic Nuclei (AGN) had their 
central (bulge) metallicities disregarded prior to determina- 
tions of the gradients. The median metallicity gradients for 
each sub-sample are given below with the number of galaxies 
in this sample following in brackets. 

• Undisturbed: -0.5017(6) 

• Disturbed: -0.4086(9) 

• Extreme: -0.2447(7) 

It should be noted here that the final figure changes dra- 
matically (to -0.5742) with the inclusion of NGC3627, which 
falls into our 'extreme' sample but shows an extremely, and 
probably unphysically, steep metallicity gradient (-2.880) 



6 Throughout this section gradients are presented in units of 
*A l °R^R2r) ' wnere ^25 is the B-band isophotal radius at a surface 
brightness of 25 mag arcsec - 2 . 
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according to our measurements. It is not known why this 
galaxy displays such a steep gradient. However, it may also 
be interesting to note that were we to include the bulge 
spectrum (initially excluded as the galaxy is classified as 
a LINER) the gradient is reduced to -2.203, leading to an 
average within the 'extreme' sample of -0.4895. 

Our data overall offer suppo rt to the much more de- 
tailed study of lKewlev et alj ([2010) , who find that the metal- 
licity gradients become shallower, or indeed disappear within 
disturbed or interacting gala xies. 



In accordance with lAnderson et al.l (|2010T ) it is also pos- 
sible in some cases to gather data on metallicity measure- 
ments at the SNe sites. Within our 'disturbed' sample there 
are 30 central CCSNe of all types (8 II+IIP and 22 SE-SNe). 
Of these we possess host HII m etallicity measurements for 
13 (from lAnderson etafl (|2010h and the recent IDS data 
presented here), all of which are classified as SNIb (7), SNIc 
(5) or SNIb/c (1). Table 3 presents the average metallic- 
ity for ea ch SNe type av a ilable here compared to r ecent 
studies bvl Anderson et all (|201Ch . lMo"diaz et all (|201ll ) and 
iLeloudas et alj ( |20 1 lT ) - The table shows that our central SNe 
have metallicities that are completely consistent with those 
presented in other studies. The SNIc within the sample are 
entirely consistent with all of the studies, w hereas the SNIb 
averag e metallicity is more consistent with ILeloudas et al.l 
|201lj ). 

The resulting conclusion from this discussion is that 
metallicity cannot drive the excess of SE-SNe within the 
central regions of interacting systems. Therefore alternate 
explanations must be pursued. 



5.2 Binary Progenitors 

The progenitors of CCSNe are wide ly discussed to arise from 
both binary and rotating stars (see lEldridge. Izzard fc Tout] 
|200S| ; iMevnet fc Maederl 120051 . respectively). The relative 
contributions from each mecha nism have been discussed ex- 
tensively in the literatur e (e.g. iPodsiadlowski. Joss, fc Hsu 
ll992j ; [Georgv et al.| [2012). Even for models where binarity is 
important, metallicity and mass m ay play an significant role 
llYoori, Wooslev fc Langeil l20ld ; lEldridge, Langer fc Tout] 
l201ll ; ISmith et alj|201lh 

However, if we are to attempt to explain the central 
excess of SE-SNe in disturbed galaxies in terms of binarity, 
we must analyse whether binary fraction can be increased 
in such environments. 

It is possible to analyse previous high density star form- 
ing regions, by looking at present day globular clusters, 
which are tho ught to be remnants of prev i ous starburst 
episo d es (e.g. Elmegreen fc Efremovl 1 19971 ; iMeurer et al.l 
19951; iLi. Mac Low fc Klessenl 120041). A recent study by 



Kruiissen et alj ( 20121 ) found that within merging galaxies 



the dynamical heating of star clusters was an order of mag- 
nitude higher in interacting galaxies than isolated ones, due 
to tidal shocks driven by the increased density. This dynam- 
ical heating is sufficient to destroy star clusters at a higher 
rate than new clusters are formed so that the total num- 
ber of stellar clusters in a merger remnant is ~ 2-50 per 
cent of the amount in the progenitor discs. However, within 
these shock heated regions the massive clusters we would 
expect to host CCSNe, are more likely to survive. Studies 
have shown that there is in fact an anti-correlation between 



binary fraction s of stars and the absolute luminosity (mass) 
of the cluster jMilone et al.ll2008l ; ISo"ilima et al.ll2007l ). 

Massive star clusters are known to have on av- 
erage higher densities an d larger velocity dispersions 
jDiorgovski fc Mevlanlll993l) . The fraction of surviving bi- 
naries is tho ught to be dom inated by binary ionisation and 
evaporation (Sollima 2008), and so within high density re- 
gions binaries are more likely to be disrupted through close 
encounters, which are both more frequen t and have higher 
mean kinetic energies (|Sollima et alJ feoiO). However, during 
encounters between binary systems and a massive single star 
the secondary star is usually ejected, therefore in clusters 
with high collisional rates there will be an inc rease in the 
number of equal mass binaries. ISollimal (|2008l ) found that 
this increased fraction is not noticeable, with an exchange 
process of only a few per cent over the entire cluster. 

Therefore this suggests that the massive clusters which 
survive galaxy interactions contain fewer binary systems. If 
these same arguments apply in the high-mass stellar regime 
then the excess of central SE-SNe cannot be explained in 
terms of an increased binary fraction. 



5.3 Modified Initial Mass Function 

In HAJ10, we characterised the modified IMF required to 
accommodate our SN observations in terms of a (radically) 
modified power-law index. However, there is to our knowl- 
edge no theoretical motivation for an IMF that is flat or even 
rising with mass, and simply changing the slope does not 
naturally explain the almost complete absence of SNIIP we 
find. Our results are more naturally accommodated by sup- 
pressing low mass SF, up to some limit that might vary from 
starburst to starburst, while leaving the form of the high- 
mass end of the IMF unchanged. To satisfy our observations, 
the requirement is that this limit should be sufficiently high 
to prevent formation of SNII progenitors in central regions 
of most (but not all) of our 'disturbed' galaxies. Above this 
limit, the IMF can take the standard form; this is uncon- 
strained by our observations. 

This concurs with the theoretical models of 
iKlessen. Spaans. fc Jappsenl (|2007l ) who find that in star- 
burst galaxies the Jeans mass increases, which affects the 
turnover of the IMF by pushing it to higher initial masses. 
This 'starburst' IMF is close to what is required to give the 
relative numbers of SE- and I IP SNe found in our study. 
The I MF calculated for the IKlessen. Spaans. fc Jappsenl 
(|2007h ' starburst' simulation would lead to a strongly 
suppressed number of SNII relative to SE-SNe, at least 
approximately consistent with the results presented here, 
under the assumption that progenitor mass differentiates 
these SN types. 

Further discussion on this form of modified IMF and its 
implications in the wider field of astronomy will be presented 
in a future paper (James et al. in preparation). 



6 CONCLUSIONS AND IMPLICATIONS 

This paper has presented the results of a re-analysis of 
HAJ10, including increased statistics and a deeper analy- 
sis of the characteristics of disturbance. The results found 
can be summarised as follows: 
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Figure 15. Metallicity gradients for a sample of the galaxies in INT IDS sample. 'Undisturbed' galaxies are in the first 
column, 'disturbed' systems in the middle column, and a selection of the 'extreme' sub-sample in the right-hand column. 



Table 3. Comparing the avera g e metallicities f or cen t ral SNIb, SNIc and SN Ibc from this study, with the literature mean 
values of lAnderson et al.l fculd ); [Modiaz et alj fculll ); iLeloudas et al.l fcOllT ). 



Study 


SNIb metallicity 


(7 


N 


SNIc metallicity 


(7 


N 


SNIbc metallicity 


a 


N 


This paper 


8.581 


0.041 


7 


8.637 


0.067 


5 


8.617 


0.035 


13 


Anderson et al. (2010) 


8.616 


0.040 


10 


8.626 


0.039 


14 


8.635 


0.026 


27 


Modiaz et al. (2011) 


8.49 


0.012 


13 


8.66 


0.010 


14 








Lcloudas et al. (2011) 


8.52 


0.05 


14 


8.60 


0.08 
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• We find a remarkable excess of SE-SNe within the cen- 
tral regions of 'disturbed' galaxies. This confirms the main 
result from HAJ10, and is an absolute excess even with re- 
spect to the star formation as traced by Hcv emission. 

• The central excess is enhanced when we compare the 
SNIbc against SNIIP populations, as expected if a sequence 
of increased envelope-stripping exists from SNIIP-IIL-IIb- 
Ib-Ic. 

• We have found no preferential 'loss' of SNIIP within the 
sample, both in terms of absolute and apparent discovery 
magnitudes, or with redshift within such a local sample of 
host galaxies. Therefore selection effects based around the 
SNe population studied here cannot drive the results found 
within this analysis. 

• The overall differences in galaxy morphology imply in- 
teraction timescales far in excess of those required to select 
the 'stripped-envelope dominated' starburst phase within 
our sample. Therefore we conclude that there are no se- 
lection effects present within our host galaxies which could 



drive the central excess of SE-SNe in disturbed galaxies pre- 
sented here. 

• Metallicity gradients within host galaxies cannot ex- 
plain this result, as our disturbed sample are likely to have 
either a small gradient, or none at all. Furthermore, the 
metallicities of the local environments of the SNIbc within 
this sample are consistent with results published data. 

• The effects of stellar rotation and binarity have been 
discussed and no evidence has been found to indicate that 
either mechanism would be increased sufficiently within our 
host environments to explain the excess of SE-SNe found. 

• Our preferred explanation is that within the central 
regions of disturbed galaxies a modified IMF exists which 
preferentially produces the most massive stars, seen here in 
through the excess of SE-SNe. 

• The radial distributions of SNIb and SNIc within 
'undisturbed' galaxies are statistically very different 
(P=0.004). This could be driven by metallicity gradients 
in these undisturbed galaxies, or radial variations in other 
properties (binarity or stellar rotation) driving envelope loss 
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in progenitor stars. This should be investigated in future 
studies. 

The results of this investigation further constrain the 
progenitors of different CCSNe subtypes, and highlight the 
need to gather as much information as possible on the spe- 
cific host galaxy environment, where direct detections are 
impossible. We maintain that in different environments the 
main mass-loss mechanism in place to produce stripped- 
envelope supernovae can differ; most notably a metallicity 
driven wind in 'undisturbed' hosts and a modified IMF in 
the central regions of 'disturbed' host galaxies. Further dis- 
cussion of the implications of a modified IMF in the central 
regions of 'disturbed' host galaxies is to follow (James et al. 
in prep). 
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